Mobility and asymmetry effects in one-dimensional rock-paper-scissors games 
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As the behavior of a system composed of cyclically competing species is strongly influenced by the 
presence of fluctuations, it is of interest to study cyclic dominance in low dimensions where these 
effects are the most prominent. We here discuss rock-paper-scissors games on a one-dimensional 
lattice where the interaction rates and the mobility can be species dependent. Allowing only single 
site occupation, we realize mobility by exchanging individuals of different species. When the inter- 
action and swapping rates are symmetric, a strongly enhanced swapping rate yields an increased 
mixing of the species, leading to a mean-field like coexistence even in one-dimensional systems. This 
coexistence is transient when the rates are asymmetric, and eventually only one species will survive. 
Interestingly, in our spatial games the dominating species can differ from the species that would 
dominate in the corresponding nonspatial model. We identify different regimes in the parameter 
space and construct the corresponding dynamical phase diagram. 

PACS numbers: 87.23.Cc,02.50.Ey,05.40.-a,87.10.Mn 



I. INTRODUCTION 

Evolutionary game theory Q, has led in the past to 
novel and interesting insights into the complex behav- 
ior that can emerge in a multispecies ecological system 
(see for example [J-Q and references therein). Cyclic 
dominance of competing "species" is not restricted to bi- 
ological systems, but can also be found in economic and 
social systems. This competition can lead to a large vari- 
ety of phenomena, ranging from the emergence of regular 
spatio-temporal pattern to chaotic dynamics. Intrinsic 
fluctuations and omnipresent nonlinearities yield a rich 
variety of scenarios, the two extreme cases being ever- 
lasting biodiversity or rapid dominance of a single species 
and extinction of the others. 

The three species rock-paper-scissors game [6T- fl9l l23j is 
one of the simplest cyclic dominance model where the in- 
fluence of fluctuations can be studied systematically. On 
the level of the mean-field rate equations this system is 
characterized by a reactive fixed point, corresponding to 
species coexistence where every species yields one third 
of the total population. When adding fluctuations, one 
has to distinguish between spatial and nonspatial sys- 
tems. In a nonspatial system, the presence of stochas- 
tic fluctuations yields as the final state one of the three 
absorbing states where only one species survives whereas 
the other two have become extinct [23| ■ Interestingly, 
the 'weakest' species, i.e. the less efficient predator, sur- 
vives in the case of asymmetric interactions [IE [24[ . In a 
spatial system without diffusion, species extinction takes 
place in one dimension [6-8] , whereas in two dimensions 



spatio-temporal pattern can emerge in the form of rotat- 
ing spirals @, EE EE EE EE El ■ Recent studies of two- 
dimensional lattice models have revealed that the mobil- 
ity of individuals, realized in the form of an exchange of 
individuals, destabilizes the spirals for larger swapping 
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rates, thus leading to species extinction in the long-time 
limit [IE El El- 
Real- world examples discussed in terms of rock-paper- 
scissors models range from coral reef invertebrates [3] to 
lizard populations in California [2(| and from competing 
bacterial strains [TJ E3 to self-organizing Min proteins 
[2l| . to name but a few examples. It is worth noting that 
some of the recent theoretical studies of low-dimensional 
rock-paper-scissors games have been motivated by the 
intriguing bacteria experiments. 

It follows from our brief discussion of the recent work 
on systems with cyclic dominance that stochastic fluc- 
tuations have important, and sometimes surprising, ef- 
fects on this type of systems. Motivated by this obser- 
vation, we wish to fully elucidate the role of fluctuations 
in low-dimensional systems composed of cyclically dom- 
inating species. We therefore propose in the following 
to complement this line of research by studying numeri- 
cally one-dimensional systems of mobile individuals, with 
species dependent interaction and swapping rates. As we 
shall show, biodiversity is made possible in one dimension 
through high mobility, which is different to the behavior 
in two dimensions where an increased mobility leads to 
species extinction. In addition, and in contrast to the 
nonspatial case, we discover for asymmetric interaction 
and swapping rates that the law of the weakest is not 
strict in the spatial case and that it depends on the sys- 
tem parameters whether the weakest survives or dies out. 

Another motivation for our study of one-dimensional 
systems comes from the recent studies of bacterial pop- 
ulations in nanofabricated landscapes (see for example 
[22} where linear arrays of coupled microscale patches 
of habitat are discussed) . These investigations open the 
intriguing possibility of a future study of competing bac- 
terial populations in quasi one-dimensional systems. 

Our paper is organized in the following way. After 
having introduced our model in Section II, we discuss in 
Section III the case of symmetric, i.e. species indepen- 
dent, rates and show that a high mobility can lead in 
one dimension to the emergence of species coexistence. 
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In Section IV we allow for asymmetric, i.e. species de- 
pendent, interaction and swapping rates and show that 
this asymmetry leads to an interesting dynamical phase 
diagram not observed in the nonspatial case. Finally, we 
conclude in Section V. 



II. ONE-DIMENSIONAL 
ROCK-PAPER-SCISSORS GAMES 

In our spatial rock-paper-scissors games we consider 
three species living on a one-dimensional support that 
are competing in a cyclic way. Calling A, B, and C 
the three species, these interactions can be cast in the 



following reaction scheme: 

A + B, B + A A + A (1) 

B + C, C + B ^ B + B (2) 

C + A, A + C^C + C (3) 



where we allow for species dependent interaction rates 
k a b, k oc , k ca . In the following we restrict ourselves to 
the case where every site is occupied by exactly one in- 
dividual, such that interactions only take place between 
nearest neighbors. In this setting mobility is realized by 
swapping the positions of individuals on neighboring sites 

mm 

A + B B + A (4) 

B + C S C + B (5) 

C + a'£ A + C (6) 

where the swapping rates s a b, Sbc, s ca can again be 
species dependent. Note that mobility can also be re- 
alized through simple diffusion steps if empty sites are 
allowed and/or more than one individual can occupy a 
given lattice site fA HI. l25j . 

In our numerical simulations, every lattice site is ini- 
tially occupied by any of the three species with the 
same probability 1/3. The total number of individuals, 
Na + Nb + Nc = N, is a conserved quantity. Here Nx is 
the number of individuals of species X in the system and 
N is the number of lattice sites. We consider periodic 
boundary conditions and use sequential dynamics. Hav- 
ing selected a neighboring pair of individuals belonging 
to species A and B, we allow for the following events: A 
dominates B with probability k ab , A and B exchange po- 
sitions with probability s Q 6, or no interaction takes place 
with probability 1 — fc a & — s Q {,. After having selected iV 
such pairs the time is increased by one unit. 

First studies of one-dimensional systems with M com- 
peting species have been published some time ago 0- 
la. [2^| . If the individuals are immobile, exact results can 
be obtained For example, it is found that the indi- 

viduals organize in single-species domains whose average 



size (A) ~ t a grows algebraically with time. The ex- 
ponent a governing this domain growth depends on the 
number of competing species as well as on the chosen 
dynamics (parallel or sequential dynamics). For the case 
most relevant to our study (M = 3 and sequential up- 
dates) the exact result is a = 3/4. Interestingly, for par- 
allel dynamics an asymptotic equivalence was established 
between the diffusion-reaction approach and the lattice 
model with immobile individuals 7] . Recent simulations 
of lattice models with sequential dynamics, where diffu- 
sion is made possible through multiple occupancy of a 
site and/or the presence of empty sites [25|, also indicate 
that inclusion of simple diffusion does not lead to a qual- 
itative different behavior. However, as we show in the 
following, this is no longer the case when the mobility of 
individuals is realized through exchanges. 



III. SYMMETRIC INTERACTION AND 
SWAPPING RATES 

Let us start our discussion of mobility effects in one- 
dimensional rock-paper-scissors games by first looking at 
the symmetric case where all rates are species indepen- 
dent. We therefore set in the following k = k ao = k bc = 
k ca and s = s ab = s bc = s ca , with k + s = 1. 

As already mentioned, earlier studies @-H, HE] of 
immobile individuals have highlighted the existence of 
species segregation through the formation of single- 
species domains that grow algebraically in time. If one 
allows for mobile individuals, the conclusions drawn in 
@, 0] remain unchanged for not too high swapping rates 
s. As an example we show in the left panel of Fig. [1] 
the space-time diagram for a system composed of 2500 
individuals where the rates are k — 0.9 and s = 0.1. Seg- 
regation sets in immediately, yielding the formation of 
numerous small single-species domains. These domains 
then coarsen, and the system rapidly decomposes into a 
few large domains. At this stage, swapping has become 
irrelevant and the behavior of the system is determined 
by the movement of the interfaces between different clus- 
ters. The time evolution of the system is then similar to 
the time evolution of the system composed of immobile 
individuals, yielding the same exponent a = 3/4 for the 
average domain size, see Fig. O Eventually, this process 
leads to the extinction of two of the three species as one 
species dominates the others in the long-time limit. 

If one keeps increasing the mobility of the individ- 
uals, one eventually reaches a critical swapping rate 
s c = 0.84(1) above which the exchange mechanism very 
effectively mixes the different species. This yields com- 
pletely different space-time plots, see the right panel of 
Fig. [TJ with a mosaic of ordered domains that emerge at 
some time and vanish again later. In addition to these 
regular pattern one also observes more chaotic space-time 
regions where no large domains are formed. This in- 
triguing spatial-temporal behavior also shows up in the 
average domain size, see Fig. [51 as (A) does no longer 
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FIG. 1: Space-time plots of one-dimensional rock-paper- 
scissors games with symmetric interaction rates k and swap- 
ping rates s. White, gray, and black correspond to the three 
different species. In the left panel, the individuals have a low 
mobility, with s = 0.1 and k = 0.9, yielding segregation of 
species, similar to what is observed for immobile individuals. 
In the right panel, high mobility, with s = 0.9 and k = 0.1, 
continuously mixes the different species, such giving raise to 
coexistence even in one dimension. The system size is 2500 
and the first 10000 time steps are shown. 



increase algebraically with time for s > s c , but satu- 
rates at a finite, s dependent, value. In this regime the 
particle densities of all three species display an irregular 
oscillatory behavior, see Fig. [3l Obviously, the mixing 
induced by the high exchange rates promotes the coex- 
istence of species, and therefore biodiversity, in our one- 
dimensional system. 

This observation of coexistence in an ecological system 
of highly mobile individuals agrees with the intuitive pic- 
ture that a mixing of species should lead to a mean-field 
like behavior. This is different to the more complex situ- 
ation in two space dimensions where for medium values 
of the mobility an intermediate regime exists in which 
rotating spirals are destabilized, yielding rapidly a uni- 




FIG. 2: Average domain size as a function of time for different 
swapping rates s, with k + s — 1. The results shown here have 
been obtained after averaging over 40000 independent runs 
for systems composed of 160000 individuals. The dashed line 
indicates the exactly known result that for s = the time- 
dependent average domain size increases algebraically with an 
exponent 3/4 for large times. 

form state where only one species survives. A coexistence 
regime similar to that observed by us in one space dimen- 
sion eventually e merg es in two dimensions for very high 
mobilities QjlEl Ef. 

Remarkably, a comparison of our results with the 
few existing results with simple diffusion (made possi- 
ble through the presence of empty sites and / or sites with 
multiple occupancy) @, HH indicates that the mechanism 
through which mobility is realized matters. At this stage 
we can only speculate on the origin of these unexpected 
differences. We note that in a field-theoretical descrip- 
tion our exchange mechanism yields a nonlinear term, in 
contrast to simple diffusion which enters linearly [27| • It 
is, however, beyond the scope of this article to explore 
this further using field-theoretical techniques. 



IV. ASYMMETRIC INTERACTION AND 
SWAPPING RATES 

There are usually no compelling reasons that require 
all interaction rates to be the same for a system with 
competing species. On the contrary, one should ex- 
pect to have species dependent rates in any real sys- 
tem, as has also been observed in a recent study of the 
self-organization of Min proteins 21]. Theoretical stud- 
ies usually neglect this aspect and exclusively focus on 
species independent interaction rates. One exception 
of which we are aware of is the study of small asym- 
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FIG. 3: Time-dependent particle densities for all three species 
for a typical run with k = 0.1 and s = 0.9 and TV = 10000 lat- 
tice sites. In this regime, where s > s c , an irregular oscillatory 
behavior is observed. 



metries in a four-state rock-paper-scissors game (where 
the forth state means that a site is unoccupied) in 
two dimensions [l8j . For that model, it was shown 
through first-order perturbation theory that rotating spi- 
rals are robust against small asymmetries. Similar con- 
clusions where drawn recently from a four-state rock- 
paper-scissors model that was constructed to describe the 
symbiosis between an ant-plant and two protecting ant 
symbionts [28]. In nonspatial rock- paper-scissors games 
an asymmetry in the reaction rates was shown to lead 
to a 'law of the weakest', as the less efficient predator 
always survives (23j . 

In principle, we have to explore for our three species 
model a six-dimensional parameter space. We will in 
the following not try to present a comprehensive study 
of all possible cases. Instead we fix the interaction and 
swapping rates between (1) B and C and (2) C and A 
particles and systematically change the rates for the pro- 
cesses involving A and B particles. Concretely, we set 
ktc = k ca = 0.4 and Sbc = s ca = 0.4. In fact, this two- 
dimensional slice through our six-dimensional parameter 
space already captures our main result, namely that the 
'law of the weakest' is not a strict one in spatial games. 

We first remark that an asymmetry in the rates yields 
the dominance of a single species both in the spatial and 
in the nonspatial game. Whereas in nonspatial games the 
weakest predator, i.e. the species X for which the inter- 
action rate k xy on another species Y is smaller than the 
other interaction rates, always survives, this is different 
in one space dimension. We discuss in Fig. 2] the survival 




FIG. 4: (Color online) (a) Survival probability of the different 
species as a function of s a b for k a b — 0.45. The total popula- 
tion is 5000 individuals. The data result from averaging over 
16000 different runs with different realizations of the noise. 
The dashed lines indicate the locations of the two transitions, 
see main text. The errors are comparable to the size of the 
symbols, (b) Survival probability as a function of time for 
k a b = 0.45 and s a t = 0.10. The different curves correspond 
to different system sizes TV, with TV = 2500 (full black lines), 
TV = 5000 (grey [green] lines), and TV = 10000 (dashed lines), 
(c) The same as (b), but now for k a b = 0.45 and s a b = 0.50. 



probabilities for the fixed interaction rate k a b = 0.45 > 
kbci k ca . Writing the survival probability of species A 
in a total population of TV individuals as Pa (TV), one 
readily verifies that the nonspatial game yields for this 
situation P B (N) — > 1 and P C (N),P A {N) — > in the 
macroscopic limit TV — > oo, i.e. the species B, which is 
dominated by A and dominates C, survives. In the one- 
dimensional case, however, the situation is much more 
complex, as the surviving species changes as a function 
of the value of the exchange rate s a b- As shown in Fig. 

for a system with 5000 individuals, we can identify 
three different regimes. Whereas for small exchange rates 
s a b < 8\ with si — 0.225(5) the species A has the highest 
survival probability, this changes for larger values as then 
either species B (for si < s a b < s 2 with s 2 = 0.32(1)) 
or species C (for s a b > S2) prevails. The transition at 
s a b = si is a smooth one, as none of the survival proba- 
bilities change dramatically when crossing the transition 
point. This is different for the transition at s a b = S2 
which is characterized by a very abrupt change of the 
survival probabilities. The latter transition has there- 
fore typical characteristics of a discontinuous transition, 
whereas the behavior at the former transition is remi- 
niscent of that at a continuous transition. We carefully 
checked the robustness of our results against a change of 
the number of individuals, see Fig. 3Jd,c. 

Repeating this for other values of fc Q b, we obtain the 
dynamical phase diagram shown in Fig. [5] Note that 
with our definition of the rates only the triangle defined 
by k a b + s a b < 1 is accessible. We can distinguish three 
different phases which differ by the dominating species: 
in phase I species A dominates, in phase II species B 
dominates, and in phase III species C dominates. These 
phases are separated by transition lines. The lower tran- 
sition line (dashed line) is characterized by the fact that 
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FIG. 5: (Color online) Dynamical phase diagram for a 
two-dimensional slice through the six-dimensional parameter 
space, with kt, c — k ca = 0.4 and sj c = s ca — 0.4. Three differ- 
ent phases characterized by the dominating species (species 
A in I, B in II, and C in III) are identified. Along the dashed 
resp. full line, one has the following relation between the sur- 
vival probabilities: Pa = Pb > Pc resp. Pa = Pb = Pc- 
The vertical dotted line indicates the line along which the 
data shown in Fig. [4^ have been obtained. The inset shows 
the transition lines for different system sizes. 



Pa = Pb > Pc, whereas along the line separating phases 
II and III we have Pa = Pb = Pc- As shown in the inset, 
the latter line does not show any size dependence. This 
is different for the line separating the regimes I and II as 
here finite-size effects are observed for the smaller system 
sizes. In fact, these finite-size dependences are compati- 
ble with a continuous transition along the line separating 
regimes I and II and a discontinuous transition along the 
line separating regimes II and III. 

Fig. [5]seems to indicate the existence at fc a & = 0.355(5) 
and s a b — 0.47(1) of a triple point where the three phases 
meet. Based on the continuous and discontinuous char- 
acter of the different transition lines, this triple point 
should in fact be a tricritical point. However, the exist- 
ing data do not allow us to fully characterize the nature of 
this triple point, due to subtle finite-size effects. Systems 
much larger than those studied here might be needed in 
order to reliable determine the location and the character 
of the triple point. 

This complex dynamical phase diagram indicates that 
the law of the weakest observed in a nonspatial game is 
not strictly valid any more for spatial systems. Thus, 
for k ao < 0.2 we indeed observe that the most inefficient 
predator eventually survives, as expected for the law of 
the weakest. However, the situation is much more com- 
plex for k a b > 0.2, as in addition to the regime where the 



weakest survives one also encounters regimes where the 
weakest dies out. 



V. CONCLUSION 

Cyclic dominance in ecological or chemical systems is 
known to yield a very rich behavior and to lead to such 
intriguing features as rotating spirals. Very recently, a 
range of new insights have been gained through the study 
of three species rock-paper-scissors games, and this in a 
spatial as well as in a nonspatial setting [l4T - [i~9l . l23j . 

In this manuscript we complement this recent re- 
search activity by a numerical study of one-dimensional 
rock-paper-scissors games with mobile individuals. We 
thereby realize the mobility through the exchange of 
neighboring individuals. In addition, we allow for species 
dependent interaction and exchange rates, thereby bring- 
ing our model closer to real systems. 

Interestingly, our study indicates that both the mobil- 
ity and the rate asymmetry yield intriguing new features. 

In presence of symmetric rates, a high mobility leads 
to a well-mixed system characterized by the prevailance 
of coexistence in one dimension. The effect of the mo- 
bility is therefore different from the effect in two dimen- 
sions where moderate mobility dissolves existing space- 
time pattern and ultimately leads to species extinction, 
whereas h igh mobility leads to the emergence of coexis- 
tence [TH, O, , similar to what we observe. Another 
surprising result is the observation that mobility realized 
through the exchange of individuals has a much stronger 
im pac t than mobility realized through simple diffusion 

For asymmetric, i.e. species dependent, rates we ob- 
serve that the spatial support induces new phenomena 
not observed in nonspatial games. Thus the law of the 
weakest, which states that the weakest predator sur- 
vives and which is a strict one in nonspatial games [23j j . 
is replaced by a dynamic phase diagram with different 
regimes that are characterized by the surviving species. 
A behavior in accordance with the law of the weakest can 
still be observed for certain ranges of the system param- 
eters, but it is no longer a strict one. 

We are not aware of any (quasi) one-dimensional ex- 
perimental realization of the rock-paper-scissors games. 
A possible way for a future realization of such a system 
could involve nanofabricated landscapes as those used 
in [22[ for a study of bacterial populations. Our study 
points out the intriguing phenomena showing up in low 
dimensions and we hope that these aspects will be ad- 
dressed in future experiments. 

In conclusion, our numerical study highlights the im- 
portance of species specific interaction and mobility rates 
in ecological systems. Our results raise some new ques- 
tions (the most intriguing being the observation that the 
effect of the mobility depends on how it is realized and 
the possible existence of tricritical points in the dynami- 
cal phase diagram) and we plan to focus on these issues 
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